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Optimal power flow (OPF) is a non-linear solution which is significantly 
important in order to analyze the power system operation. The use of 
optimization algorithm is essential in order to solve OPF problems. 
The emergence of machine learning presents further techniques which 
capable to solve the non-linear problem. The performance and the key 
aspects which enhances the effectiveness of these optimization techniques 
are compared within several metaheuristic search techniques. This includes 
the operation of particle swarm optimization (PSO) algorithm, firefly 
algorithm (FA), artificial bee colony (ABC) algorithm, ant colony 
optimization (ACO) algorithm and differential evolution (DE) algorithm. 
This paper reviews on the key elements that need to be considered when 
selecting metaheuristic techniques to solve OPF problem in power 
system operation. 
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1. INTRODUCTION 

In recent decades, concentrated research effort were employed to solve complex optimization 
problems. Several methods are integrated to ensure the techniques able to operate in optimal rate. One of it is 
optimal power flow (OPF) problem in power system operation. The major concern is to reduce total 
operation cost while satisfying various constraints in the system. Computational intelligence assists to 
understand the problems with various perspective by introducing several metaheuristic techniques which are 
inspired by biological systems [1]. 

Evolutionary computation is used to find optimal solutions. This includes genetic algorithm (GA) 
which is used for search and optimization purposes; evolutionary algorithm (EA) which is excellent on 
decision making; and multi-agent based algorithm which works best on sorting and usually used for data 
mining. Traditionally, Newton-Raphson method that is used to identify the optimal values which is started by 
making an initial guess. This method requires large number of iterations to achieve optimal solution for the 
problems. Modern metaheuristic algorithm improves the process by using randomization to improve the 
convergence rate. Convergence of an algorithm is crucial to achieve optimal results which might take long 
time to converge. While reducing the time will lead to premature converge which means the solution 
provided does not guarantee global optima for the problem that being considered [2]. 

Several population-based heuristic search techniques are used to provide high quality optimal 
solution in acceptable timescale. Several changes in parameter setting are made to suit the nature of the 
optimization technique such as particle swarm optimization (PSO) algorithm, ant colony optimization (ACO) 
algorithm, artificial bee colony (ABC) algorithm, firefly algorithm (FA) and differential evolution (DE) 
algorithm. These algorithms used the iterative process and self-organize system to search the space of 
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solution using evolutionary operation. Since it is impossible to create a universal algorithm which able to 
solve variety of optimization problems, researchers have come out with several metaheuristics with similar 
framework that contains several changes in the parameter. Although it is good to have variety of nature 
inspired algorithms such as glowworm swarm optimization (GSO) algorithm, monkey algorithm (MA), 
lion optimization algorithm (LOA), grey wolf optimization (GWO) algorithm and etc, the main focus is the 
needed to achieve a balanced exploration and exploitation process to achieve complete convergence with 
reasonable time based on the complexity of the problems. This paper reviews on the key elements that need 
to be considered to solve OPF problem in power system operation based on five population-based 
metaheuristic search techniques which are PSO, FA, ABC, ACO and DE. 


2. OPTIMAL POWER FLOW PROBLEM FORMULATION 

The combination of power flow equation and economic dispatch equation can be simplified into 
non-linear function. Standard OPF problem can be formulated to minimize the objective in the system and 
satisfy system equality and inequality constrains as: 


Minimize f (x,u) 
Subject to g (xu) <0 

h (x,u) = 0 (1) 
where : x : The dependent variables 

u : The independent variables 


f(x,u) : Objective functions 
g (xu): Set of inequality constraints 
h (x,u)_: Set of equality constraints 


2.1. The Dependent Variable (State Vector) 
x’ = [Po, Vir ... Vint, Qai ... Qane, Su... Sint] (2) 


where: PGI : Generator active power at slack bus 
VL : Load-bus voltage 
QG : Generator reactive power 
S1 : Transmission line flow 
NL : Number of load buses 
NG : Number of generators 
nl : Number of transmission lines 
T : Transformer tap setting 


2.2. The Indipendent Variable (Control Vector) 


ul = [Pa2...Pona, Vai... VGna, Qc1.--Qene, T\... Tr] (3) 
where: PG : Generator active power at PV buses (except at the slack bus PG1) 

VG : Generator voltage at PV buses 

QC : Shunt VAR compensation 

NT : Number of regulating transformers 

NC : Number of VAR compensators 


2.3. The Equality Constrains (Typical Load Flow Equation) 


NB 
Po — Pr, — V; »y [Gy cos(6, — 5, ) 
j=i 


+ B,sin(d; — 6 )| = 0 (4) 
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NB 


Qai — Qi - Vi x V; [Gy cos(6; — 8 ) 
j= 


+ B,sin(S, — 6 )| = 0 


(5) 
where: NB : Number of buses 
PD : Active load demands 
QD : Reactive load demands 
Gi : Transfer conductance 
Bij : Susceptance 
2.4. The Inequality Constrains (Operation Limits of Equipment) Generator constraints: 
Vai" < Vai < Vai ,i=1, ...,.NG (6) 
pwn < Po < pmax ,i1=1,...,NG (7) 
QGi" S Qa S VWI ,i=1, ..., NG (8) 
Transformer constraints: 
je < T; < Tj" x : i= 1, a NT (9) 
Shunt VAR constraints: 
OB < Qci $ QM ,i=1, ....NC i 
Security constraints: 
Vi" s Vis s VE, i=1,...,NL ay 
’ yMaX 5 
$5 i=1,..., nl i 


The OPF problems solution are designed by considering both normal and contingency operation 
conditions. This is to ensure that contingencies factor will not cause any violations of limits. Since the 
inequality constrains are self-limiting therefore, the optimization techniques using evolutionary computation 
techniques can be performed to overcome the problem. 


3. PERSPECTIVE OF OPTIMIZATION ALGORITHM 
3.1. Exploitation (Intensification) and Exploration (Diversification) 

Exploitation or intensification in an algorithm is important to generate new solutions which better 
than existing solutions referring to the problem of interest. This allow high convergence rates but usually 
trapped in local search. However, the local search is limiting the end solution since it did not compare the 
information globally and it depends on the starting location. Exploration or diversification on the other hand, 
search for information globally and able to generate the solution with enough diversity. Since it needs to 
search in large space, this will slow down the convergence rate [1]. 

The balance between exploitation and exploration will be able to converge the solution faster and 
the potential to achieve global optima is high. The current algorithms use different parameter setting which 
cause certain algorithm to achieve fast convergence rate or able to achieve global optima. Randomization 
also been used by certain algorithms to make sure the convergence rate is fast. Since there are so many 
algorithms that perform very well in very specific problem, No Free Lunch (NFL) theorems is introduced. 
NFL states there is no universal algorithm that efficient enough for every problem. This is because the main 
purpose of OPF is to obtain high quality solution within acceptable timescale. So there is no reason to obtain 
average solution for all problems that being considered. However, some researchers suggest that NFL did 
exist on some types of problems such as coevolutionary approaches [3]. 
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3.2. Evolutionary Operators 

GA and EA uses genetic operators such as recombination (crossover), mutation and selection 
(survival of fitness) [4]. This allows diversity from mutation and crossover process; and also from the 
selection and survival of fitness high quality information can be obtained. This is important to solve complex 
and discontinuous problems. The rate of crossover and mutation will influence the global search rate 
(exploration) while the selection increase the exploitation action. In general, crossover is used to mix within a 
subspace and allowing to reach convergence of system. Mutation is a randomization technique which makes 
the iterative process for local search faster [2]. The selection is used to allow the system to evolve towards 
desired states. This happen by allowing the current solution compared with previous best solution, if the 
current solution generate better solution, the previous solution will be discarded. This repetitive process is 
carried out until the final best solution is achieved. Taxonomy of Evolutionary Computation (Optimization 
Algorithm) is shown in Figure 1. 


Evolutionary “ 
o-ea:- 


Iterative Method 
Gradient 
Newton's Method 


Nearest Neighbor Heuristic y / 


Insertion Heuristic : " 
rd Genetic Algorithm 
, / Evolution Strat 
Soest / Differential Evolution 
Threshold Algorithm Cc P 
/ Evolutionary Programming 
Biogeography-based Optimizer 
Learning Classifier System 
Neighborhood Search Method — Scatter Search 


Enumerative Method 





Saving: Heuristic 





Particle Swarm Optimization 
Ant Colony Optimization Algorithm 
Firefly Algorithm 
Artificial Bee Colony Algorithm 
Termite Algorithm 
Wasp Swarm Algorithm 
Monkey Algorithm 
Glowworm Swares Optimization Algorithm 
. Lion Optimization Algorithm 
~ Grey Wolf Optimization Algorithm 
reer mmen Bat Aigoridkee 
Hunting Search 
Krill Herd Algorithm 
Sayre Echolocation 
yhale Optimization Algorithm 






Teaching Learning Based Optimization 
Harmony Search 
Tabu Search 
Group Search Optimization 
Imperialist Competitive Algorithm 





Simulated Annealing 
Gravitational Local search 
Big-Bang Big-Crunch 
Gravitational Search Algorithm 
Changed System Search 
Central Force Optimization 


a Weed Optimization Algorithm 
Artificial Immune System 
Bacterial Foraging Optimization 
Dendritic Cell Algorithm 


Interior Search Algorithm 
Mine Blast Algorithm 
Soccer League Competition Algorithm 


Artificial Chemical Reaction Optimization Algorithm 
Black Hole Algorithm 
Ray Optimization Algorithm 
Small World Optimization Algorithm 
Galaxy Based Search Algorithm 
Curved Space Optimization 





Exchange Market Algorithm 
Group Counseling Optimization Algorithm 


Figure 1. Taxonomy of Evolutionary Computation (Optimization Algorithm) 


3.3. Heuristic Techniques 

Heuristic technique is an approach to solve problem or learning using practical method to obtain 
satisfactory solution within reasonable time. This occurs by using similar rules from data set which has been 
trained to solve similar problems. Based on Figure 1, heuristic techniques can be classified into 4 categories 
which are heuristics, metaheuristics, matheuristics and hyperheuristics. 

Matheuristics is a hybrid techniques that combines exact and approximate method to solve 
mathematical problems of interest using their exploitation features. These techniques uses improved 
mathematic programming techniques based on metaheuristics and vice versa [5]. This technique also shows 
significant improvement in local searching and local branching during neighborhood exploration [6], [7]. 

The approximate techniques which involves heuristics, metaheuristics and hyperheuristics applies 
exploitation and exploration. However, for heuristics method, it focus more on the exploitation which allows 
the problem able to obtain good solution in reasonable time. But they are usually trapped in local optima due 
to minimal exploration. 

Metaheuristics techniques are suitable to optimize large and complex data. This is possible due to 
the behavior which accepts a temporary deterioration of the solution. Active exploration and exploitation 
allows it to randomly explore to obtain global optimum solution. However, the true balance of exploration 
and exploitation in any algorithms or problems still does not have a clear line and is supported with NFL 
theorem. 
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Hyperheuristics focus more on exploration compared to exploitation. The set space also not as usual 
space set of solution. This is because hyperheuristics is used as heuristics to search for optimal heuristics to 
be used. It is used to handle several classes of problem rather than solving just one problem [8]. The idea is to 
compensate the strength and weakness of known heuristics to select suitable method to be used for the 
assigned problems. 


3.4. Classification of Metaheuristics 

There are four main criteria that are used to classify various technique of metaheuristics. Firstly, the 
method used for every iteration is classified between trajectory and discontinuous methods. Usually 
trajectory methods is used for single point search in which the process just closed walk on the neighborhood 
graph. Discontinuous methods commonly used by population based techniques which allows it to make large 
jumps in the neighborhood graph for high exploration in the search space. 

Second criteria is the involvement of memory in the algorithm. Since most of the techniques which 
is based upon either short term memory or long term memory, hence there is no clear group for this criteria. 
The purpose of having short term memory is to avoid cycling in similar solution by forbid the revisiting to 
the recent solution. Long term memory on the other hands, allows good diversification and intensification. 

Third criteria is the way the algorithm interact with neighborhood structures. It is either kick-moves 
operation which interact with single neighbor or using operation that able to interact with various 
neighborhood such as mutation and crossover. Final criteria that been consider is the objective function 
during operation of the technique. Some algorithm uses static objective function to achieve optimal solution. 

However, there is also algorithms which uses reward or penalty to modify the objective function. 
This dynamic objective function allows the search process to be guided in the correct direction to obtain 
global optimal solution. All four criteria are used to determine the main trait for each algorithm then 
classified into the group as shown in Figure 1. 


4. POPULATION-BASED METAHEURISTIC SEARCH TECHNIQUES 

The basic operation and theory to obtain OPF have been discussed. This section will focus on the 
exploration and exploitation. Techniques on current popular algorithms such as Particle Swarm Optimization 
(PSO), Firefly Algorithm (FA), Artificial Bee Colony Algorithm (ABC), Ant Colony Optimization (ACO) 
and Differential Evolution Algorithm (DE). Those five algorithms are compared due to the versatility in 
optimization process in various problem and the development in various study and hybridization techniques. 


4.1. Particle Swarm Optimization Algorithm (PSO) 

PSO was developed by Kennedy and Eberhart in 1995 which is based on the behavior of bird and 
fish schooling in nature [9], [10]. PSO is widely used due to the details explanation in algorithm operation in 
their book ‘Swarm Intelligent’ [11]. This allows PSO to be easily hybridized with other algorithms to 
compensate with the PSO limitation when dealing with several problems. The movement of particles depend 
on stochastic component and deterministic component. This two main components did used to determine the 
current global best location and current best location. Since the movement is random, active comparison 
occurred and allowing global best solution can be obtained over time. The vector can be represented as: 


vit? = vj + aeylg* — xi] + Beglxj - xj] (13) 
vitt = xf + vit? (14) 
where: Xi : Position of a particle 
Vi : Velocity of a particle 
1 €2 : Random vectors that carry values 
between 0 and 1 
a,p : Learning parameters or 
acceleration constants, typically 
azxp 2 
g : Current global best solution 
Xj : Individual best 


Based on (13) and (14), PSO consists of mutation and selection as whole. Crossover is not present in 
PSO algorithm. This allows a high exploration rate in order for fast convergence rate. Since no crossover 
happens during the process, it is high possible that the current global best solution has prematurely 
converged. 
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However, there is possibility to achieve global optimal solution. The randomization for initial 
location allowing for the process take different route. Multiple random location will possibility lead to global 
optimal solution. 


4.2. Firefly Algorithm (FA) 

FA was developed by Xing-She Yang in 2008 which based on the flashing patterns and behavior of 
tropical fireflies [12]-[14]. The brightness of the flashing will be used for better exploration. The vector can 
be represented as: 


xj") = xf + Boe Vi (xj — xi) + ae; 


(5) 


where: i : Movement of a firefly 
J : More attractive (brighter) firefly 
Bo : Attractiveness at r=0 
a : Randomization parameter 
€f : Vector of random numbers drawn 


Based on (15), FA is easier to be implemented and flexible. The mutation that used from random 
numbers drawn from the Gaussian distribution and Levy flight allowing local and global search. The rate of 
searching can be controlled with randomization parameter. When compared with PSO, FA did not have 
global best solution which means FA did not undergo selection process. 

However, FA uses ranking system which will be updated during the loops. The attraction feature is 
unique for FA which allows it to automatically classifies the solution into several groups [15]. This due to 
fact that the local attraction is stronger than the attraction from another location. The classification allows FA 
suitable to be used as multi objective problems and can be applied in various applications which are highly 
non-linear problems. 


4.3. Artificial Bee Colony Algorithm (ABC) 

ABC algorithm was developed by Karaboga in 2005 which is based on social behavior of bees [16]. 
ABC algorithm is divide into three group which are forager bees (employed), observer bees (onlooker) and 
scout bees. ABC did used the concentration of solution as their indicator which known as the amount of 
nectar intake from different food sources. Scouts and employed is used for mutation process to generate 
randomization then the solution will be selected by the onlooker. Employed bees will forage every food 
sources available with ratio of one bee per food source. When the food sources area been discard, forager bee 
will become scout to search for new food sources area by referring to the information given by onlooker 
bees. When compared to honey bee algorithm, ABC has advantage when the forager bees converted to scout 
directly. The vector for forager bees phase, observer bees phase and scout bee phase can be represented 
respectively as: 


vy = Xi + Diy (Xiy — Xj) (16) 
fit 
y= san oc 
7 SN fi 
fit 
sc (17) 
—_ min max min - 
xy = xi" + (xf xj™") rand(0,1) (18) 
where: vj; : Source position 
Sy : Size of employed or onlooker bees 
x; : D-dimensional vector (i =1, 2,..., Sy) 
. : Number of optimization parameters 
(/€1,2,...,D) 
- : Comparison between two food 
7 positions 


k : Random chosen indexes different 
from i (k € 1, 2,..., Sy) 

Oi; : Random number between -1 and 1 

xj": lower bounds of parameter j 

xj"“* > upper bounds of parameter j 
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Based on (16)-(18), there is no crossover to maximize their randomization. This allowing the 
algorithm to adapt with tough optimization since the operation is simplify with the specific task which 
already been assigned for the particular roles [17], [18]. ABC in general have the ability to globally search 
which cause the convergence rate relatively slow and the exploitation ability is limited due to the absence of 
crossover process. 


4.4. Ant Colony Optimization Algortihm (ACO) 

ACO algorithm was developed by Dorigo in 1999 which based on social behavior of ant [19-21]. 
When compared to FA which using attraction, ACO used similar techniques as ABC which uses the 
concentration of route (quality of solutions). This concentration shows the quality of the solution which is 
generated for discrete combinatorial problems. The vector can be represented as: 


a 4b 
Dime ij ij 
Uo on a 4B 
Lj=1 Oi dij 


(19) 


where: i : Probability in network at particular 
node (before) 
J : Probability in network at particular 
node (after) 
i; : Concentration on the route between i 
and j 


d : Desirability of the same route 


Based on (19), ACO did used evaporation over time features to avoid premature convergence. 
The main technique by ACO is the probability of choosing a route and the evaporation rate [22]. The best 
solution will display with high concentration of route due to the number of route used is higher due to 
traveling time is short, thus less amount of evaporated within that short period. This short term memory 
allowing the best route remain higher in concentration without causing any redundant from another local 
maximum which been discovered previously. 


uy 


4.5. Differential Evolution Algorithm (DE) 

DE algorithm developed by Rainer and Kenneth in 1997 which based on creating new off-spring 
from parent chromosomes instead of classical crossover or mutation [23]-[25]. DE did used crossover, 
mutation and selection. The unique of DE is the self-adaptive mechanism and selection process cause by their 
mutation scheme. DE used greedy selection process that provide slightly advantages of converging 
performance for the new and parent solution that win the selection [24], [26], [27]. DE can be stated into four 
stages which are initialization, mutation, crossover and selection. 

During the initialization, the component can be represented as: 


x, (0) = x} + rand(0,1) (x/! - x}) 60) 


where: rand(0,1)  : Uniformly distributed random 


number 
xj' : parameter of upper bound 
x} : parameter of lower bound 


i : Population members 

J : Component of population 

While the mutation phase, donor vector is created from the changes of population member which 
can be represent as: 


vit +1) = Xp (OO + F(X 2) — X30) 
(21) 
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where: Xi) : Changes of population member 
Vict) : Donor vector 
Xrl, : three distinct points taken randomly 


X2,X3 from the current population and not 
coinciding with the current xi 
F : Scalar number which between 0.4 
and 1.0 
The crossover process did have two crossover schemes that commonly used which are exponential 
and binomial. This allowing diversity of population. Next, the process that ensures the population size 
remains constant is known as selection stage. The selection process is determined by the concept of survival 
of fitness which gives advantage for the new solution and their parents. This process can be represented as: 


x(t+1)= u,(t) if (fF (u(t) s f (%;,(t)) 


x(t+1) = x(t) if (f (x(t) < f (u(t) 
(22) 


where u, j(t) : The crossover information 
fQO : Function to be minimized 


5. CONCLUSION 

The main key for evolutionary operators to ensure optimal result are crossover, mutation and 
selection. Crossover is efficient in exploitation. While with the extensive usage of mutation and selection 
enables the ability of exploration. The current global best helps to speed up the convergence but there is 
probability for premature convergence. Therefore, the parameter setting for the algorithm is crucial to ensure 
the balance between exploitation and exploration which will improve the overall performance of the 
algorithm. The complexity of the algorithm also needs to be considered as a form of user friendliness. 
The research on the parameter setting for various problem is still insufficient. This can be done by 
understanding the NFL theorem on problems that it is applied on. The hybridization of several algorithms is 
useful to enhance the rate and result of the convergence. 
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